Spin-i fermions on spin-dependent optical lattices 
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We study the phase diagram of one dimensional spin-i fermionic cold atoms. The two "spin" 
species can have different hopping or mass. The phase diagram at equal densities of the species is 
found to be very rich, containing Mott insulators and superfluids. We also briefly discuss coupling 
ID systems together, and some experimental signatures of these phases. In particular, we compute 
the spin structure factor for small momentum, which should allow the spin gap to be detected. 
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Quantum engineering Q, 0] of strongly correlated 
many-body systems has recently become possible thanks 
to the spectacular advances in trapping ultracold atoms 
in optical lattices 0, S IE @ or in microchip traps • 
This has led to the study of models that would otherwise 
be hard to realize in solids, which may shed light on ba- 
sic issues in quantum many body physics, including the 
understanding of e.g. the origin of high-T c superconduc- 
tivity in doped copper oxides. In particular, correlated 
boson HHE@, Bose-Fermi [Hj, and Fermi [HE El 
systems have received much experimental and theoretical 
attention in recent times. 

In sharp contrast to electrons in solids, in cold atomic 
systems, different types of atoms (different hypcrfmc 
states or different atomic species) can be trapped and 
controlled independently, such that the hopping, strength 
and sign of interactions (inter- or intra-species) and den- 
sities can be continuously tuned. For example, Mandel 
et al. H managed to control independently the periodic 
potential for each atom type loaded in an optical lattice. 
A p-wave Feshbach resonance 01 can create a tunable 
asymmetry in the interactions in a multi-species Fermi 
gas. All this, of course, leads to a much richer physics, 
which remains to be understood. 

In this paper, motivated by these recent developments 
and the availability (now or soon) of fermions in elon- 
gated traps la, |jfl , we study the interesting effects of hav- 
ing different Fermi velocities for two species of fermions 
in one dimension (ID). With equal densities of the two 
species, this system is different from the case of a two- 
leg spinless ladder Jljl. l2l| or the spin-i electrons in a 
magnetic field (lil l2l| . One main result of this paper 
is the phase diagram as a function of velocity difference, 
for equal densities (Fig. GJ. With repulsive interactions, 
a finite velocity difference breaks the SU(2) spin sym- 
metry and turns the gapless Tomonaga-Luttinger liquid 
(TLL) into an Ising spin-density wave with a spin gap. 
Demixing may occur if one type of fermions has a very 
tiny velocity. With attractive interactions, a singlet su- 
perfluid (SS) of bound pairs of fermions of different types 
may give way to a charge density wave (CDW) of pairs 



for sufficient velocity difference. We also briefly study 
the effects of a small tunneling term coupling an array 
of ID tubes together. In particular, if there are different 
densities of fermions in neighboring tubes, a triplet super- 
fluid (TS) may become stable for repulsive interactions. 
Finally, we describe how to detect the spin gap in these 
phases by measuring the dynamical spin structure factor, 
which we have computed. In contrast to previous stud- 
ies j^, El , we have worked out the phase diagram for 
equal number of spin up and down fermions as a function 
of the (Fermi) velocity difference and considered coupling 
the ID systems together. 

We study the following generalized Hubbard model: 



H = - t„ (4 m c CTm+i + H.c.) +u^2>, 
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(1) 



This Hamiltonian describes a ID Fermi gas with contact 
interactions (related to U) prepared with Nq fermions 
(i.e. we work in the canonical emsemble) and loaded 
in a fD optical lattice with m = 1, . . . , M lattice sites; 
n am = c^mCam, and a =|, 4 is the spin index that may 
refer to two hyperfine states, or to e.g. 6 Li and 40 K. 
Even though there may be no true spin symmetry, we 
will continue to use the spin (and magnetic) language 
to describe this binary mixture. We assume the number 
of fermions of each spin species is separately conserved, 
i.e. one spin type cannot be converted to another. Moti- 
vated by the experimental considerations above, we allow 
for different hopping t a for different spins. This Hamil- 
tonian may be realized in cither a quasi- fD chip trap 
|7| or in a 2D optical lattice, which is made up of an 
array of f D gas tubes 0, Q weakly coupled by a hop- 
ping t± <C min{i|,ij}. When the finite "charging en- 
ergy" (due to the finite length of each tube) exceeds the 
renormalized hopping Ej cx t± 0, Il7l | , the tubes are de- 
coupled from one another and a set of independent ID 
tubes is recovered ||j. Although we assume there is a 
(spin-dependent) periodic potential parallel to the tubes 
such that applies, much of what is discussed below 
also applies in the absence of this potential when the two 
species have different masses. More discussion on engi- 
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FIG. 1: Schematic phase diagram for the model in Eq. Q 
with equal number of spin up and down fermions away from 
half-filling (i.e. 7V T0 = N l0 / M/2) . The interaction 
strength is U and a = \tf— tj,|/(tf All phases (SDW: spin 
density wave, CDW: charge density wave, SS singlet super- 
fluid, TS triplet superfluid) exhibit a spin gap A s (however, 
A s = for U = and z = with (/ > 0. A cartoon of the 
type of order characterizing each phase is also shown. In the 
area between dashed lines the dominant order (either CDW 
or SS) depends on the lattice filling (see text for more detailed 
explanations). In the SG phase, spin up and down fermions 
are segregated (demixed). 



neering Hamiltonians like 0} can be found in 

We first study the homogeneous ID system in the ther- 
modynamic limit. Finite-size and trap effects will be dis- 
cussed below in connection with possible experiments. 

The weak coupling limit \U\ <C min{£j,^} can be 
solved by taking the continuum limit of (JTJ in the stan- 
dard way |2lj | and linearize the dispersion around the 
Fermi points ±k F = irNo a /Ma, (M is the number of 
lattice sites in the tube). This leads to the so-called "g- 
ology" representation |2l| with a finite number of cou- 
pling constants representing low energy scattering pro- 
cesses. The coupling gZ< (<72_i_) is the scattering ampli- 
tude for processes where a small momentum q is ex- 
changed between fermions of equal (opposite) spin at 
opposite Fermi points, for arbitrary values of k F . On 
the other hand, gi± is the back-scattering amplitude 
where two fermions of opposite spin exchange a momen- 



vp — ^">p 

./Vof = Nq^; <?3_l is the amplitude for umklapp scattering 



turn q w 2k p — 2k F = 2k p, and is relevant only when 



(q w 2k F + 2k F = tt /a) and is important only at half- 
filling: iVoj + N^i = M . Thus for generic fillings, qi i 
and g3_L are irrelevant, and the system is a TLL |2l|), 
which has a completely gapless spectrum of two distinct 



branches of phonons. 

We focus here on the case iVoj = iVoj. ^ M/2. The 
case of a half filled lattice iVof = iVoj. = M/2 is more 
involved and will be reported elsewhere [lj| 

The physical properties can be established by analyz- 
ing the renormalization group (RG) flow of the various 
scattering amplitudes upon the varying of a cutoff such 
as the temperature T . To second order in the interaction 
parameters, the RG flow is [l7|: 



2/2|| 



-crJ/lJ_ , 



V2± 



r lV l 2\\ 



2?/2_L Vl± 



(2) 
(3) 



where y a = g a /Trhv are dimensionless couplings, v = 
(f| + V[)/2 the mean velocity and r a — v/2v a \ y a = 
dy a /d£, with £ = ln(A/T). Eqs. l|2l:ill can be mapped to 
the RG equations of the Berezinskii-Kosterlitz- Thouless 
(BKT) transition in terms of y2 S = — Yl^ r o-?/2|| 2?/2_l 
and The behavior of the BKT equations is entirely 
determined [2lJ by the constant of motion C 



yL/2(r T n + 1) = (Ua/hv) 2 z 2 /{2 - z 2 ), where z = |t T - 
£jj/(tf is the key velocity difference parameter. For 
z = we recover the well-known results for the spin- 
symmetric Hubbard model |2lj |. However, for z ^ and 
U ^ 0, C > 0, the scattering amplitude yi±(£) diverges 
as the system is cooled down to its ground state. This 
signals the formation of bound states, and the opening 
of a gap in the spin sector (the charge excitations remain 
gapless). For z 1, the gap has thus the characteristic 
BKT form A, - A e~ A l^c ^ A e — ^4'/l*r-*t. I , where A is 
of the order of £f rj t± and A, A' are constants. Note that 
this gap is non-perturbative in |tj — t±\. 

The properties of the spin-gapped phase depend on the 
sign of U. Ground states of ID systems are characterized 
by the dominant form of order that they exhibit, which 
is typically quasi-long range in character, true long-range 
order being only possible in ID when a discrete symme- 
try is broken. For U > and z ^ 0, then yi±(£) —* oo, 
and a bosonization study shows that the dominant 
order is a spin-density wave (SDW) and the subdominant 
order is triplet superfluidity (TS). In the attractive case 
(U < 0), as z is increased, the dominant order changes 
from a singlet superfluid phase (SS) to a charge density 
wave (CDW), with CDW and SS being the subdominant 
order in the former and latter case, respectively. We 
wish to point out that our analysis takes fully into ac- 
count the marginal (in the RG sense) coupling between 
the gapless charge and the gapped spin modes arising at 
z ^ 0, which leads to an often substantial decrease in 
the value of the Luttinger-liquid parameter K c (propor- 
tional to the charge compressibility). In particular, for 
U < we find that K c goes from K c > 1 to K c < 1 as 
z is increased, which changes the character of the domi- 
nant correlations from SS to CDW, as described above. 
A summary of the phase diagram is shown in Fig. ^ 
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The weak-coupling regime smoothly crosses over to the 
strong coupling regime \U\ 3> max{t|, tjj, as confirmed 
by a strong coupling expansion of ■ We only give here 
the main steps, technical details can be found in [l7| . 
It is simplest to first consider a half-filled lattice with 
7V|o = N~io = M/2. For the strongly repulsive case U » 
max{f| , } fermions cannot hop around and there is a 
gap of order U to charge excitations. Degenerate pertur- 
bation theory then shows that in this limit the Hamil- 
tonian in (JJJ maps to the Heisenberg-Ising (XXZ) spin 
chain H X xz = JJ2 m [ S ™ ' s m+i + 7 S'm^+i] , where 
the S m denotes the spin operator at site m, J — it^ti/U, 
and the anisotropy 7 = (if — tj) 2 /2tjij cx z 2 . Thus, 
for unequal hopping (z > 0), the chain is in the Neel 
phase (SDW with true long-range order), and has a 
spin gap which for small z is A s ~ J e - " 4 A/7 = 
J e~ A /l*T _t 4-l [2^. Note the same non-perturbative de- 
pendence on t| — ij as for the weak coupling regime. 
Away from half-filling, the system is described by a t-J- 
like model with anisotropic spin interactions. The charge 
gap is destroyed (K c = \ close to half-filling H^]), but the 
spin gap remains and the dominant order is still SDW. 
Physically, the finite velocity difference breaks the SU(2) 
spin symmetry to the lower Z2X U(l). Thus, the TLL 
becomes an Ising anti-ferromagnet in the spin sector. 

For U <C 0, degenerate perturbation theory shows that 
P|l is equivalent to a model of tightly bound fermion 
pairs (hard-core bosons) annihilated by b m — C| m C| m . 
Their hopping amplitude is J = 4t^ti/\U\, and they in- 
teract with strength V = 2(t| + tJ)/\U\ when sitting at 
nearest-neighbor sites. This model can be mapped to the 
above XXZ chain via b m — > S~ and (b^bm — |) — > S* n . 
At half-filling, charge excitations are gapless for equal 
hopping and SS is the dominant order 21]. However, 
with unequal hopping the spectrum of the tube is fully 
gapped, becoming a CDW with true long-range order, 
a spin gap of order \U\ (energy to break a pair), and a 
charge gap A c ~ J e _yl /l iT_tj -L Away from half-filling, 
the spin gap remains ~ \U\ but the bosons are able to 
hop (i.e. the charge gap disappears). It is worth pointing 
out that very close to half-filling for z ^ 0, the dominant 
order is CDW since K c — > ~ , as can be inferred from the 
exact solution of the XXZ chain [2ll |25|. However, as 
the filling deviates more and more from half-filling, K c 
rises above one and the system becomes a ID superfluid 
(SS). This change in the character of the dominant order 
also takes place at constant filling, provided the system 
is sufficiently far from half-filling: a SS (K c > 1) can 
turn into a CDW (K c < 1) as \z\ is varied at strong cou- 
pling. This agrees with the above weak coupling analysis. 
Note that at very low density (N 0cr /M — > 0), and at least 
for not too different velocities, only a SS phase is possi- 
ble: in this limit, Q maps to a continuum (Gaudin-Yang 
like) model of interacting fermions with spin-dependent 
mass. For U — > —00, the fermions pair up to become a 
ID superfluid (SS) of tightly bound pairs with irrelevant 



residual interactions between the pairs. 

Finally both the weak and strong coupling analysis de- 
scribed above break down for sufficiently large |tj — tjj. 
For weak coupling, linearization of the free fermion dis- 
persion is no longer justified when tf t± (or vicev- 
ersa), that is, for z — > 1. In the large \U\ limit degener- 
ate perturbation theory becomes quite subtle. Unfortu- 
nately, rigorous results are available only for t-\ = or 
t± = (z = 1), which is the limit of the Falicov-Kimball 
model. In ID, Lemberger 0] (see also |2jj) has proved 
that spin up fermions segregates from spin down ones for 
U > U c > at equal densities. There is no se greg ation 
for U < at equal densities. As argued in \VA l2d| it is 
quite likely that this segregated phase will survive also 
when \z\ is not one but close to one. 

The predicted phase diagram of Fig. 2] for a single ID 
tube can be directly tested experimentally in cold atoms. 
However, it is also interesting to analyze the case when 
the tubes are weakly coupled by tunneling between the 
tubes. We thus briefly describe the phase diagram for 
an array of coupled ID tubes in a 2D optical lattice ge- 
ometry 0, The Hamiltonian for each tube at site 
R of the 2D lattice is as in Eq. (JTJ, with all fermion 
operators now carrying the R label. The hopping be- 
tween the nearest neighbor tubes at R and R' is de- 
scribed by H± = -t_L E<R,R') Em^cIlte^R'mi where 
tj_ -C min{i|,t|}, but such that fermions can now over- 
come the 'charging energy' of the finite-size tubes. In 
general, when the isolated tube has a gap A s <C t±, H± 
is a relevant perturbation (in the RG sense) that leads to 
coherent hopping of fermions from tube to tube. Thus 
the ground state will most likely be a very anisotropic 
3D Fermi liquid, which in turn may become unstable to 
3D CDW/SDW formation or 3D BCS superfluidity under 
appropriate conditions. This limit has been much studied 
in the past e.g. in connection with organic superconduc- 
tors (see [24| for a review) . We shall not consider it here, 
and instead we study if ^ so that the gap A s 3> t± . 

Since the tubes (or at least a large number of them near 
the center of the trap due to inhomogeneity effects) can 
develop a sufficiently large spin gap as described above, 
coherent hopping between tubes is now suppressed. 
However, the term H± can generate, through virtual 
transitions which are second order in t± , intcrtube inter- 
actions of two kinds |2l|: i) particle-hole pair hopping 
generates spin-spin and density-density interactions: 
Hi = X/m,(R,R') [«Al Sr to • Sr'„, + Vj_n m R,n m R/], 
and ii) fermion pair hopping yields H% = 
J c±E m ,(R,R') b R m b Km+j, where V±,JcL ~ t 2 JA s , 
b-Rm — c |R m c lRm and ja < £ s a distance smaller than 
the spin correlation length, £ s cx Aj 1 . The dominant 
term then drives a phase transition to a 3D ordered 
phase 17]: for the U < case, if the tubes are in 
the SS phase, then the dominant process is fermion 
pair tunneling, and the tubes develop 3D long-range 
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superfluid order. The low-temperature properties of 
this system become identical to the superfluid of bosons 
studied in However, if the tubes are in the CDW 
(U < and sufficiently large z) or in the SDW (U > 0) 
phases, the dominant interactions arise from hopping of 
particle-hole pairs and lead to insulating phases that are 
either 3D CDW or SDW. The ordering temperatures in 
all cases (at small t±) are power-laws: T c oc A s (ij_/A S ) Q , 
with oT 1 = 2(2 — d) and d the scaling dimension of 
the dominant inter-tube interaction. Interestingly, the 
SDW or CDW ordering is anisotropic: incommensurate 
(relative to the optical lattice) along the tube, but 
commensurate perpendicular to the tube direction. 

Particle-hole hopping may drive a transition to a 3D in- 
sulating state with density wave order only if the density 
in neighboring tubes are equal or very similar: for a par- 
ticle and a hole to hop coherently at low-temperatures, 
they must be extracted from opposite Fermi points of one 
tube and must match the momenta in the neighboring 
tube by momentum conservation. If the mismatch in the 
density between tubes is sufficiently large, particle-hole 
hopping is suppressed and only the hopping of fermion 
pairs (which carry zero net momentum) is possible. The 
system will then order as a superfluid. Interestingly, for 
U > 0, TS is the subdominant order in the spin gapped 
phase of the tubes. Suppression of particle-hole pair hop- 
ping may then lead to a 3D triplet superfluid. We note 
that T c for these cases is also a power law of t±/A s . 

The phase diagram shown in Fig. ^ holds, strictly 
speaking, in the thermodynamic limit. In real 2D op- 
tical lattices only a finite number of fermions (~ 10 5 |6() 
can be loaded, but we expect all predicted phases to ap- 
pear. Due to the finite size, the phase boundaries will 
not correspond to true phase transitions, but rather to 
sharp crossovers. The trap can lead to phase coexistence 
and even to suppression of quantum criticality j2(|, but 
we are concerned here with the phases themselves and 
not with the quantum critical points between them. 

The most important signature of the single tube phases 
that we predict is the existence of a spin gap, A s . By 
measuring the absorption of a laser that causes Raman 
transitions between the two hyperfine states a =T,1, it 
should be possible to measure [28j the small momentum 
limit of the dynamic structure factor S s (q, w), which is 
the Fourier transform of S s {r, r', t) = (S+jr, t)S~(r', 0]). 
Using the so-called form factor approach |29( , we find [l7| 
that at T < A s , the structure factor rises from zero as 
y / (haj) 2 - (2A S ) 2 for Hlu > 2A S . 

Concerning the coupled tubes, the most exotic phase 
is of course the triplet superfluid (TS). To 'engineer' it, 
we need to suppress particle-hole hopping by making the 
number of fermions in neighboring tubes sufficiently dif- 
ferent. This could be achieved by imposing a rapid spa- 
tial variation of the trap potential, or better, by means 
of a biperiodic optical potential in the direction perpen- 
dicular to the tubes. The coherence properties of the 



3D superfluid phases could be probed by exciting low 
frequency collective modes in the transverse direction to 
the tubes. Coherent oscillations should exist only in the 
superfluid phases and not in the insulators. 
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